Finite thickness of shear bands in frictional viscoplasticity and implications for lithosphere dynamics

Permanent deformations in the lithosphere can occur in the brittle as well as in the ductile domain. For this reason, the inclusion of viscous creep and frictional plastic deformation is essential for geodynamic models. However, most currently available models of frictional plasticity are rate independent and therefore do not incorporate an internal length scale, which is an indispensible element for imposing a finite width of localized shear zones. Therefore, in computations of localization, either analytical or numerical, resulting shear zone widths tend to zero. In numerical computations, this manifests itself in a severe mesh sensitivity. Moreover, convergence of the global iterative procedure to solve the nonlinear processes is adversely affected, which negatively affects the reliability and the quality of predictions. The viscosity that is inherent in deformation processes in the lithosphere can, in principle, remedy this mesh sensitivity. However, elasto‐viscoplastic models that are commonly used in geodynamics assume a series arrangement of rheological elements (Maxwell‐type approach), which does not introduce an internal length scale. Here, we confirm that a different rheological arrangement that puts a damper in parallel to the plastic slider (Kelvin‐type approach) introduces an internal length scale. As a result, pressure and strain and strain rate profiles across the shear bands converge to finite values upon decreasing the grid spacing. We demonstrate that this holds for nonassociated plasticity with constant frictional properties and with material softening with respect to cohesion. Finally, the introduction of Kelvin‐type viscoplasticity also significantly improves the global convergence of nonlinear solvers

Post-orogenic extension and metamorphic core complexes in a heterogeneous crust: the role of crustal layering inherited from collision. Application to the Cyclades (Aegean domain)

International audienceThe development of metamorphic core complexes (MCC) corresponds to a mode of lithospheric continental stretching that follows collision. In most of the models that explain the formation of the MCC, high thermal gradients are necessary to weaken the lower crust and to induce its ascent. Such models fail to explain the exhumation of high pressure-low temperature metamorphic rocks in metamorphic core complex structures as observed in the Cycladic Blueschists in the Aegean domain. Besides, account for the lithological crustal stratification induced from collision has never been tested. In this paper, we use fully coupled thermomechanical modelling to investigate the impact of structural heritage and initial thermal gradient on the behaviour of the post-orogenic continental lithosphere. The models are designed and validated by petrological, structural and time data from the Cyclades. As a result, high thermal gradients (Moho temperature higher than 800°C) are neither necessary nor always sufficient to induce the development of a metamorphic core complex. At the contrary, the rheological layering of the crust inherited from collision is a first-order parameter controlling the development of extensional structures in post-orogenic settings. 'Cold' MCC can develop if the crust is made of a strong nappe thrust on top of weaker metamorphic cover and basement units, as observed in the Cyclades

Toward robust and predictive geodynamic modeling : the way forward in frictional plasticity

Strain localization is a fundamental characteristic of plate tectonics. The resulting deformation structures shape the margins of continents and the internal structure of tectonic plates. To model the occurrence of faulting, geodynamic models generally rely on frictional plasticity. Frictional plasticity is normally embedded in visco‐plastic (V‐P) or visco‐elasto‐plastic (V‐E‐P) rheologies. This poses some fundamental issues, such as the difficulty, or often inability, to obtain a converged equilibrium state and a severe grid sensitivity. Here, we study shear banding at crustal‐scale using a visco‐elasto‐viscoplastic (V‐E‐VP) model. We show that this rheology allows to accurately satisfy equilibrium, leads to shear band patterns that converge upon mesh refinement, and preserves characteristic shear band angles. Moreover, a comparison with analytic models and laboratory data reveals that V‐E‐VP rheology captures first‐order characteristics of frictional plasticity. V‐E‐VP models thus overcomes limitations of V‐P and V‐E‐P models and appears as an attractive alternative for geodynamic modeling

Formation of metamorphic core complex in inherited wedges: A thermomechanical modelling study

International audienceMetamorphic Core Complexes (MCCs) form when a thickened domain with a low-strength lower crust is submitted to extension. These structures are characteristic of post-orogenic extension, and field observations suggest that several MCCs rework a crustal nappe-stack emplaced before extension begins. These MCCs therefore develop within heterogeneous crusts that contain pre-existing dipping heterogeneities, such as thrust faults and dipping nappes in a crustal wedge. Although very common, this first order structural inheritance has never been considered in studies modelling MCCs. Our contribution therefore investigates the effect of an inherited crustal wedge structure on the dynamics and kinematics of formation of the MCCs, using fully coupled thermomechanical modelling. The wealth of petrological, structural and time informations available in the Cycladic MCCs (Aegean domain) allows setting up more realistic initial conditions for the experiments than usual flat-lying setups. It also allows the results of the numerical computation to be directly validated with final geometries, P-T paths and exhumation rates. The experiments using dipping heterogeneities are characterised by a much more complex evolution and final structure than their flat-lying layered equivalents. Dipping heterogeneities drive lateral strength contrasts and help to re-localise the deformation on successive detachments. The dip of the inherited wedge structures imposes kinematic constraints on the flow, which provides a model that explains the regional scale asymmetry of the Cycladic MCCs. The P-T paths, the exhumation rates and the final crustal structure that come out of an initial shallow-dipping wedge model provide a much more realistic comparison with their natural counter-parts than common flat-lying models. Other parameters, like crustal-scale density inversion, thermal structure and creep law parameters are of second order when compared to the initial wedge structure. Being little dependent on these second order parameters, the proposed model for the formation of MCCs within inherited crustal wedges is likely to be applied to other areas where the MCCs formed in a nappe stack involving continental basement

Burial and exhumation in a subduction wedge : mutual constraints from thermo-mechanical modelin and natural P-T-t data (Sch. Lustrés, W. Alps)

The dynamic processes leading to synconvergent exhumation of high-pressure low-temperature (HP-LT) rocks at oceanic accretionary margins, as well as the mechanisms maintaining nearly steady state regime in most accretion prisms, remain poorly understood. The present study aims at getting better constraints on the rheology, thermal conductivity, and chemical properties of the sediments in subduction zones. To reach that goal, oceanic subduction is modeled using a forward visco-elasto-plastic thermomechanical code (PARA(O)VOZ-FLAC algorithm), and synthetic pressure-temperature-time (P-T-t) paths, predicted from numerical experiments, are compared with natural P-T-t paths. The study is focused on the well constrained Schistes Lustrés complex (SL: western Alps) which is thought to represent the fossil accretionary wedge of the Liguro-Piemontese Ocean. For convergence rates comparable to Alpine subduction rates (∼3 cm yr−1), the best-fitting results are obtained for high-viscosity, low-density wedge sediments and/or a strong lower continental crust. After a transition period of 3-5 Ma the modeled accretionary wedges reach a steady state which lasts over 20 Ma. Over that time span a significant proportion (∼35%) of sediments entering the wedge undergoes P-T conditions typical of the SL complex (∼15-20 kbar; 350-450°C) with similar P-T loops. Computed exhumation rates (<6 mm yr−1) are in agreement with observations (1-5 mm yr−1). In presence of a serpentinite layer below the oceanic crust, exhumation of oceanic material takes place at rates approaching 3 mm yr−1. In all experiments the total pressure in the accretionary wedge never deviated by more than ±10% from the lithostatic component

Strain localisation in mechanically layered rocks beneath detachment zones: insights from numerical modelling.

International audienceFully dynamic numerical simulations have been designed in order to asses how the orientation of mechanical layering in rocks controls the orientation of shear bands and the depth of penetration of strain in the footwall of detachment zones. Two parametric studies are presented. In the first one, the influence of stratification orientation on the occurrence and mode of strain localisation is tested. The second parametric study shows that results are length-scale independent and that orientation of shear bands is not sensitive to the viscosity contrast or the strain rate. Based on the results, a conceptual model for strain localisation under detachment faults is presented. In the early stages, strain localisation occurs at slow rates by viscous shear instabilities but as the layered media is exhumed, the temperature drops and the strong layers start yielding plastically, forming shear bands and localising strain at the top of the shear zone. Once strain localisation has occured, the deformation in the shear band becomes extremely penetrative but the strength cannot drop since the shear zone has a finite thickness

Granite intrusion in a metamorphic core complex: The example of the Mykonos laccolith (Cyclades, Greece)

International audienceThe Aegean domain is a well-suited place to study the formation of metamorphic core complex (MCC) and to investigate the role of syn-tectonic granites on their development. In the northern Cyclades, the Mykonos-Delos-Rhenia MCC is characterised by the intrusion of a kilometer-scale Late Miocene pluton of I-type granitoids within a migmatitic gneiss dome. New combined AMS (Anisotropy of Magnetic Susceptibility) and microstructural studies on the Mykonos granitoids together with recently published thermochronological data allow us to use the granitoids as strain markers. The Mykonos granitoids form a laccolith-like intrusion with a N70°E long axis. The laccolith is strongly asymmetric with an outlying root zone to the SW and a major body mainly developed to the NE. The laccolith construction is due to successive pulses of more or less differentiated magma that intruded the Cycladic Blueschist Unit. The attitude of stretching markers suggests an important (about 60°) vertical-axis local rotation phenomenon in the cycladic upper crust during the exhumation of the Mykonos MCC. Structural data suggest a four-stage evolution of the Mykonos MCC: (i) a first stage characterized by flat shearing toward the N-NE and by the formation of a domal structure in migmatitic paragneisses with multi-scale generation of folds with axes either perpendicular or parallel to the regional stretching, as a result of the interplay between regional N20E°-directed extension and EW shortening; (ii) a second stage marked by the emplacement of the Mykonos laccolith at 13.5 ± 0.3 Ma at the top of the migmatitic paragneisses; (iii) the third stage corresponding to the development of protomylonitic foliations and lineations in the whole laccolith in high to medium temperature conditions; (iv) the late stage marked by an acceleration of the exhumation of the Mykonos MCC. This exhumation was accommodated by important rotations of upper crustal blocks. During the end of the exhumation processes, around 10 Ma, deformation localized at the top of the laccolith in semi-ductile conditions and then in brittle conditions in the major detachment plane. Our study shows that the Cycladic plutonism event had no role on the initiation of the MCC. However, the geometry of the Mykonos intrusion supports that the magmas are "sucked" into the direction of regional extension and that the intrusion of magmas has caused an acceleration of the last stages of the MCC development. This acceleration was marked by a very fast exhumation of the laccolith after its emplacement

The North Cycladic Detachment System

International audienceLow-angle normal faults accommodate a large part of continental post-orogenic extension. Besides the intrinsic rheological characteristics of the continental crust that may lead to the formation of shallow-dipping shear zones at the brittle–ductile transition, the role of pre-existing low-angle structures such as large thrusts has been proposed by several authors. We explore this question with the example of the North Cycladic Detachment System (NCDS) that is composed of a series of distinct detachments cropping out on the islands of Andros, Tinos and Mykonos, separating the Cycladic Blueschists in the footwall from the Upper Cycladic Nappe in the hanging wall. We show that these extensional structures are part of a single large-scale structure (more than 200 km along strike) that reactivates the Vardar suture zone. It extends eastward on Ikaria and westward offshore Evia and Thessalia where it probably connects to recent shallow-dipping normal faults evidenced on published seismic reflection profiles. The NCDS started its activity in the Oligocene concommitantly with the Aegean extension, and was still active in the Late Miocene. It has exhumed a series of metamorphic domes from southern Evia to Mykonos below low-angle detachment systems, made of low-angle normal faults and low-angle ductile shear zones. The ductile shear zones and the faults were created with a low dip and they kept the same attitude throughout their exhumation. We identify three main detachments that are part of a continuum of extension on the NCDS : Tinos detachment, Livada detachment and Mykonos detachment. A fourth detachment (Vari detachment) is the reactivation of an Eocene exhumation-related structure. Deformation in the footwall is characterized by intense stretching and flattening. Using the spatial evolution observed along strike from Andros to Mykonos we construct a history of formation of the NCDS starting with the reactivation of former thrusts leading to the exhumation of high-temperature metamorphic domes. The Aegean example shows that reactivation of earlier shallow-dipping discontinuities can play a fundamental role in continental post-orogenic extension

Thermal imprint of rift-related processes in orogens as recorded in the Pyrenees

International audience19 The extent to which heat recorded in orogens reflects thermal conditions inherited from 20 previous rift-related processes is still debated and poorly documented. As a case study, we 21 examine the Mauléon basin in the north-western Pyrenees that experienced both extreme 22 crustal thinning and tectonic inversion within a period of ~30 Myrs. To constrain the time-23 temperature history of the basin in such a scenario, we provide new detrital zircon fission-24 track and (U-Th-Sm)/He thermochronology data. The role of rift-related processes in 25 subsequent collision is captured by inverse modeling of our thermochronological data, using 26 relationships between zircon (U-Th-Sm)/He ages and uranium content, combined with 27 thermo-kinematic models of a rift-orogen cycle. We show that the basin recorded significant 28 heating at about 100 Ma characterized by high geothermal gradients (~80°C/km). Our 29 thermo-kinematic modeling and geological constraints support the view that subcontinental 30 lithospheric mantle was exhumed at that time below the Mauléon basin. Such a high 31 geothermal gradient lasted 30 Myr after onset of convergence at ~83 Ma and was relaxed 32 during the collision phase from ~50 Ma. This study suggests that heat needed for ductile 33 shortening during convergence, is primarily inherited from extension rather than being only 34 related to tectonic and/or sedimentary burial. This should have strong implications on tectonic 35 reconstructions in many collision belts that resulted from inversion of hyper-extended rift 36 basins

Influence of basement heterogeneity on the architecture of low subsidence rate Paleozoic intracratonic basins (Reggane, Ahnet, Mouydir and Illizi basins, Hoggar Massif)

The Paleozoic intracratonic North African Platform is characterized by an association of arches (ridges, domes, swells, or paleo-highs) and low subsidence rate syncline basins of different wavelengths (75–620&thinsp;km). The Reggane, Ahnet, Mouydir and Illizi basins are successively delimited from east to west by the Amguid El Biod, Arak-Foum Belrem, and Azzel Matti arches. Through the analysis of new unpublished geological data (i.e., satellite images, well logs, seismic lines), the deposits associated with these arches and syncline basins exhibit thickness variations and facies changes ranging from continental to marine environments. The arches are characterized by thin amalgamated deposits with condensed and erosional surfaces, whereas the syncline basins exhibit thicker and well-preserved successions. In addition, the vertical facies succession evolves from thin Silurian to Givetian deposits into thick Upper Devonian sediments. Synsedimentary structures and major unconformities are related to several tectonic events such as the Cambrian–Ordovician extension, the Ordovician–Silurian glacial rebound, the Silurian–Devonian Caledonian extension/compression, the late Devonian extension/compression, and the Hercynian compression. Locally, deformation is characterized by near-vertical planar normal faults responsible for horst and graben structuring associated with folding during the Cambrian–Ordovician–Silurian period. These structures may have been inverted or reactivated during the Devonian (i.e., Caledonian, Mid–Late Devonian) compression and the Carboniferous (i.e., pre-Hercynian to Hercynian). Additionally, basement characterization from geological and geophysics data (aeromagnetic and gravity maps), shows an interesting age-dependent zonation of the terranes which are bounded by mega-shear zones within the arches–basins framework. The old terranes are situated under arches while the young terranes are located under the basins depocenter. This structural framework results from the accretion of Archean and Proterozoic terranes inherited from former orogeny (e.g., Pan-African orogeny 900–520&thinsp;Ma). Therefore, the sedimentary infilling pattern and the nature of deformation result from the repeated slow Paleozoic reactivation of Precambrian terranes bounded by subvertical lithospheric fault systems. Alternating periods of tectonic quiescence and low-rate subsidence acceleration associated with extension and local inversion tectonics correspond to a succession of Paleozoic geodynamic events (i.e., far-field orogenic belt, glaciation).</p